Cosmological dynamics of D-cceleration 
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We study the cosmological dynamics of an effective theory for a strongly coupled scalar field in 
the moduli space of A/" = 4 supersymmetric Yang-Mills theory recently proposed by Silverstein and 
Tong, called "D-cceleration". We discuss various Energy Conditions in this theory. Then we prove 
the inflationary attractor property using the Hamilton- Jacobi method and study the phase portrait 
as well as the cosmological evolution of the scalar field. 

I. INTRODUCTION 

^P ■ It is a great challenge to implement inflation in string/M-theory. Currently, to achieve this goal, there are mainly 

^~^ \ two types of models: one is the modification of the Friedmann equation in the braneworld scenario (see, e.g. [l| for a 

' review); another is specific model of inflaton field from string/M-theory (see, |3| for a recent important development). 

C In this paper, we will consider a scenario belonging to the second type recently proposed by Silverstein and Tong 0]. 

. ^ , In 13] , Silverstein and Tong considered the infiaton as the scalar field in the moduli space of the A/" = 4 supersym- 
metric Yang-Mills (SYM) theory. They call this model "D-cceleration". In the strong coupling region, the effective 

^T^ description of the theory is via gravity and string theory by using the AdS/CFT correspondence |a|- Then in the 

' gravity side, the causal speed-limit in the bulk will translate into a speed limit on moduli space. This has also been 

,__! \ observed before by Kabat and Lifschytz |j]. Interpreted in the field theory side of the correspondence, this result 

^ ■ reflects the breakdown of the moduli-space cr-model approximation due to the growing importance of the higher 

S^O '■ derivative terms. The dynamics will be changed dramatically from the canonical scalar field theory due to this speed 

C . limit and may lead to inflations without a flat potential. 

\l ■ The effective action comes from, as written in field theory variable, the Dirac-Born-Infeld action of a probe D3-brane 

^^ . moving in AdSsxS^ SIBl with metric signature {— , -I-, -I-, -I-}, 

^ . 9ym J ^ ^'X^ ^^ '■^' X 

1^ where gvM is the Yang-Mills coupling, and A = Nqym is the 't Hooft coupling, in which N is the rank of the SYM 
; I ■ theory. The string coupling gs is given by gs — 5y^,^/(47r). It is worth commenting that in the pure Af — 4: SYM 
"^ I theory, there does not have the potential V{(l)) (The dynamics of the 4> field without the potential is analyzed in 
Ci ' Sec. 3. 2 of PI)- The potential V{(j)) arises only after couple to other sectors involved in a full string compactification. 
IJ ' However, as we will show in Sec^ a simple argument based on Strong Energy Condition shows that if we want the 
• i-H , (j) field to be a reasonable candidate of infiaton, we must consider a non-zero potential. 

^ ' It is interesting to note that when A(9^(/))^/0'* ^ 1, expanding the action (Q) to first order in A(9^(/))^/^^, it just 

H ] reduces to the lagrangian of a canonical scalar field theory £ — ^{d^cj))'^ + V{(p). So the lagrangian Jp) can be viewed 
■ - - ■ as a correction to the canonical scalar field theory when the kinetic term (9^0)^ is large. On the other hand, the 
lagrangian (QJ is also similar to the effective action of tachyon field !q| £ = V{T)y/l + {d^J^Y, which is a physical 
example of the general scenario of "k-infiation" y\. So it is also_conceivable that the dynamics of the theory Q will 
be similar in some respects to the tachyon field. The author of |3| mentioned that the theory ^ can be viewed as a 
physical example of k-inflation. However, the general lagrangian of k-inflation is of the form C = K {(j))p{X) , where 
X — ^dfj_(t>d^(t>, K and p are two arbitrary functions. Then we can see that the lagrangian iQJ actually does not 
belong to this form. So this lagrangian is not an example of k-inflation in a precise meaning and the analysis of the 
dynamics of k-inflation cannot be applied to it. 

The most notable feature of the lagrangian (^ is that it automatically imposes a restriction on 0, the speed limit: 



^ = -3^/d^-y^[^Jl + ^^q^ + n0)-^], (1) 



A 



(2) 
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Thus when (j) is small, the (j) field will be enforced to roll slowly without the need of a flat potential. 

The energy density p and pressure p of the cj) field follows from the action ^ (Note that since we are interested in 
the cosmological dynamics of the (j> field, it is reasonable to assume that cj) is spatially homogeneous) , 

V{q^) - 5^, (3) 



a/i^aJv^ 



^ ;i-A<;6V<^4_y(<^) + r_, (4) 

and the equation of motion for cj) can be obtained either by varying in Q or by the energy conservation equation 
p + 3H{p + p) = 0, 

^ _ 5|! , i£ , 3H(1 - f )* + (V'W -f^l- p-" = (5) 

It is easy to check that 



^/A 



(6) 



t + V\/(t)i 

where <f>i is the initial value, is always a solution of the evolution equation ((SJ which satisfies the speed limit exactly, 

i.e. (/)^ = ^, irrespective of the precise form of the potential V. However, this cannot be a stable solution: from 
the expression for the energy density p, i.e. Eq.J2|), we can see that when (f) satisfies the speed limit exactly, the 
energy density tends to infinity (This also has an interesting dual interpretation from AdS/CFT: when the speed of 
the D-brane tends to unity, according to relativity, its energy tends to infinity) . This means that the action (^ is not 
sensible when (j) satisfies the speed limit exactly. 

II. ENERGY CONDITIONS 

When facing a new field theory that intended to apply to cosmology, one natural question one would like to ask is: 
under what conditions it will satisfy or violate the various Energy Conditions j^l? 

First, it is easy to see from Eqs.Q, l@J that (j) always satisfies the Weak Energy Condition: p > and p + p>i). 
This implies that the equation of state of w = p/ p always satisfies w > — 1. Thus this theory will not suffer from 
the possible instability of vacuum [2j that might cause problems in a phantom field theory which is characterizes by 
UJ < -1. 

Second, with a little more labor, it is also easy to say that the Dominant Energy Condition: p > \p\, is also always 
satisfied by 0. Thus, from the vacuum conservation theorem of Hawking and Ellis j^, (see also [lO| a recent review 
and a simplified proof), the energy of cj) field cannot propagate outside the light-cone. 

At last, let's come to the Strong Energy Condition: p + p> and p + 3p > 0. From Eqs.Q, |Q, 



20' 
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P + 3p=-^il-\Jl- A0V04) + ^ _ - 2^(0). (7) 

^ ^1 - A0V04 

Then it is easy to see that when V = 0, (j) will always satisfy the Strong Energy Condition. As is well-known, for 
a matter to drive an accelerated universe, it is just necessary and sufficient for that matter to violate the Strong 
Energy Condition. So if 1^ = 0, as in the case of pure Af = A SYM theory, (f> field cannot be a viable candidate of 
infaton. Whether cj) field can violate the Strong Energy Condition depends on the precise form of the potential. In 
this sense, the (j) field is more like a canonical scalar field rather than a tachyon field. Since in the latter case, whether 
the tachyon field can violate the Strong Energy Condition depends only on the value of |T| (see, e.g. [III). 

III. INFLATIONARY ATTRACTOR PROPERTY 

Inflation can be predictive only if the solution exhibits an attractor behavior, where the differences between solutions 
of different initial conditions rapidly vanish (see Sec. 3.7 of [TjI for a review). Thus this section is devoted to the 
study of the inflationary attractor property of the (f) field for a general potential V{4>). 



The Friedmann equation reads, 

where H = a/a is the Hubble parameter. 

We wiU use the Hamihon- Jacobi formulation for the Modified Friedmann equation [l3| , which is a powerful and very 
effective tool for analyzing the inflationary attractor property , 14| . This formulation is also very useful in obtaining 
solutions of the evolution equations |j3|- In this formulation, we will view the scalar field (j) ^ the time variable and 
this requires that the (j) field does not change sign during inflation. Without loss of generality, we can choose > 
in the following discussions. If this is not satisfied, it can be brought about by redeflning (j> -^ —(j). 

The Hamilton- Jacobi formulation of the Friedmann equation reads |3J 

1/(0) = ■i{gsMl)H{<f>f ~ ^'^ ^Jl/{glM^) + ^XH'ict^Y/r + ^ (9) 

and 



^l/(g2Af4)+4Ai/'(0)2/04 



(10) 



Supposing Ho{(f>) is any solution to Eq.®, which can be either inflationary or non-inflationary. We consider a small 
homogeneous perturbation 5H{(j)) to this solution. The attractor property will be satisfled if it becomes smaller as 
increases. Substituting H{(f)) — Ho{(j)) + SH{(t)) into Eq.Q and linearizing, we find that the perturbation obeys 



2H;,6H' = 3Ho^l/{g^sM^) + A\H'{^f /4>HH (11) 

which has had the general solution 

5H{4>) = SH{^,) expif iH,^l/{9lM^) + A\H'{4>Y /^^^] (12) 

where 5H{4>i) is the value at some initial point (j)i. From Ea. H10(l . since H'^ and d4) have the opposite sign, if Hq is an 
inflationary solution, all linear perturbations are damped at least exponentially. 

Thus we conclude that the inflationary attractor property holds for the fleld 0. See Fig.l for the phase portrait of 
the potential (|13|) where we can clearly perceive this conclusion. 

IV. PHASE PORTRAIT AND COSMOLOGICAL EVOLUTION 

In this section, we will study numerically the dynamics of the (p fleld for the simplest potential that may lead to a 
power law acceleration Q 

V{c^) = im202 (13) 

It has been shown in Q that at late time, 4> will just tend to solution ©. We choose this potential is because that 
in addition to its simplicity, as argued in 3], higher power terms in the effective potential will be small compared to 
the quadratic term without signiflcant tuning. 

However, the lagrangian Q cannot be trusted for arbitrarily small because of the back reaction of the probe 
brane on the geometry. It is shown in 3] that the theory can be trusted only for greater than a critical value [3J, 

</>2 > lg,m2 ^ cfi, (14) 

Also in p] , it has been argued that the potential (|13|) can drive a power law inflation only if m satisfies 



where the Planck scale M^ = (SttG)"^. We have verified numerically that this is indeed the case. Thus in the 
following discussions, we will impose this bound on the parameter m. 

In the following discussion, we will fix the tension parameter A — 100 throughout for illustrative purposes. 

To study the evolution, it is convenient to rewrite the evolution equation (jSJ as a set of two first order equations 
with two independent re-scaled variables Yl = (j)/M and Y2 — 0/M^ where M is the inflationary energy scale. Then 
in terms of li and I2 , the speed limit (|2Jl and the bound (|14|l translate into 

Yi < ^ (16) 



A 



and 



Y^' > \gsrl ^ Fic (17) 

where r^ = m/M < 1 for any reasonable value of the scalar field's mass m. 

We define a dimensionless parameterri = M/Mp for describing the energy scale fraction. In the following discussion, 
we will always take r^ = 10^^, which corresponds to the GUT scale for the M. Then the evolution equation can 
be written as an dynamically autonomous system: 

Yo 

Y; = —^ (18) 

^ H/M ^ ' 



>^2 = T7^[^ - ^ - 3^(1 - ^)Y, {rlY, 1^)(1 - ^f'-] (19) 



AKj 
H/M' Yl A "M'^ Yf' ^'™'' ^T^^' ~Y^ 

where a prime denotes differentiation with respect to the number of e-folding N defined by 



iV = In — = / Hdt, (20) 



thus dN = Hdt = d\na. 

And in terms of Yi and Y2 , the Fricdmann equation ||SJ) can be rewritten as 






Also in terms of r^ and r^, the bound p5|l translates into 




(22) 

In the following we start the numerical analysis of the dynamical system. Fig.l shows the trajectories in the phase 
plane (^1,12) with g^ = 1 x 10~* and r^ = 105s/(rf A) = 10"^. Yi will take values in the region Yic < Yi < 1; F2 will 
take values in the region — Yj^/VA < Y2 < Y^ jyfX. We can see that there is a curve that attracts all the trajectories 
which correspond to the slow-roll curve. This confirms our conclusion in the last section for the specific potential Ijl^ll . 
For all the allowed initial value, tends to the slow-roll curve quite soon after it begins to evolve. For simplicity, let 
us then analyze in detail the dynamics of (/> for small initial Y2 that satisfies ^^il^u ^ 1- Since Yn ^ 1 and A = 100, 
this corresponds to Y2i ^ 10~^. 

See the left panel of Fig.l where we have taken Y^i ~ 10^"^. By the numerical calculations, the dynamics of </> can 
be divided qualitatively into three stages. First, when it begins to evolve, it quickly evolves to the slow-roll curve. 
Second, it undergoes a slow-roll stage until Y\ ~ 0.2. Third, it evolves as the speed Hmit curve. 

In the first stage, from the left panel of Fig.l, we can consistently set Yi = const, and because we have assumed 
a small initial Y2, we can set 1 - XY^/Y^ « 1. Thus from Eq.(|2Tl), H/M = const, and {H/Mf ~ IQ'^Y^. Because 
Y2 ~ 10~3 and Fi - 1 in this stage, we have Yi/{HYi/M) < -Y2 with Y2 < 0. Thus Eq.(|Tnil reduces to 

Yi - -iY2 - ^n (23) 
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FIG. 1: Phase portrait for the field (/> in the re-scaled coordinates (Yi, Y2) for ps = 1 x 10~* and Vm = (10gs/{rf\)y^^. The left 
panel show the behavior for small initial value of ¥2- The right panel shows the behavior for initial value of Y2 close to the speed 
limit bound. In both panels, the dotted line corresponds to the curve that satisfies the speed limit exactly, i.e. Y2 — —Y^j^fX. 
The dashed line in the right panel corresponds to the slow-roll approximated equation 12711 . 



which can be integrated to give 



i"2 = {y2^ 



2,H/M 



Yi)cxp(-37V)- 



iH/M 



Yi 



(24) 



Thus in this region, I2 is damped exponentially to the slow-roll curve with Y\ kept almost constant as we can see that 
from the upper right of the left panel. Furthermore, we can see that during the slow-roll stage, Y^ ~ ~r'^-^/{iH/M). 
In the second stage, since \Y2 /Y^ ^ lO^'^, we may set 1 — XY2 /Y^ « 1 and Y2 « in the dynamical system as 
described by Eqs. (|18|l . ()19|l . and (|21|l . Then the system can be approximated by 



and 



6K2 H , 






This is a simple algebraic system of (Fi, I2) and can be solved exactly to give 



Y2 = [3— - \h{ — Y + 24r2 1^ 



(25) 



(26) 



(27) 



The dashed line in the left panel of Fig.l corresponds to the phase portrait of this slow-roll approximated evolution 
equation. Furthermore, since (iJ/M)^ ~ \Q~^Y^ and r,^„ — 10^'^, Ea. lf77|l can further reduce to 



Y. 



iH/M 



Yi 



(28) 



This is consistent with Ea. (|24|l . 

In the third stage, when Y\c <Y\ < 0.2, the attractor curve coincides with the limit curve that satisfies the speed 
limit exactly, i.e. Y2 — —Y^l\f\. This has already been pointed out in 3] as the late time behavior of by using the 
Hamilton- Jacobi method. However, as we have showed, will behave in this way well before it reaches the critical 
point (^c- 



CONCLUSIONS AND DISCUSSIONS 



In this paper, we studied the cosmological dynamics of the D-cceleration scenario. We proved the inflationary 
attractor property using the Hamilton- Jacobi method and studied the phase portrait. Thus in this respect, this 



model can be a reasonable model of inflation. Of course, there are still many works needed to be done such as 
confrontation with observation to see whether this is a viable model of inflation. 

It is interesting to ask whether the theory lU can also be a viable candidate for quintessence which can drive an 
acceleration in the current universe (see, e.g. (15i| for a recent review). This is interesting because the model building 
of quintessence faces the same problem as inflation: it requires the quintessence field to be in a slow-roll phase which 
requires a flat potential. But this is rather difficult to achieve from a particle physics point of view. So the important 
feature of D-cceleration, i.e. slow-roll without a flat potential, is also very attractive as a quintessence candidate. 
However, the problem is that, assuming we still work with the potential (|13|l . the theory will make sense only under 
the condition that (j) > <f>c- So if we want to apply this model to the late universe, we must check that will satisfy 
(j) > 4>c after inflation for a wide range of initial conditions. Although we still cannot present a definite answer to this 
question, through numerical computation, we found that the evolution of <p will be extremely slowed down before it 
reaches (pc- So this possibility still remains. 

Following the same line of reasoning after the effective theory of tachyon is proposed, it is also interesting to consider 
the evolution of the D-cceleration in brane- world scenarios such as the well-studied Randall-Sundrum II model ^^. 
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